Electronic properties of curved graphene sheets 
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A model is proposed to study the electronic structure of slightly curved graphene sheets with 
an arbitrary number of pentagon-heptagon pairs and Stone-Wales defects based on a cosmological 
analogy. The disorder induced by curvature produces characteristic patterns in the local density of 
states that can be observed in scanning tunnel and transmission electron microscopy. 
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The recent synthesis of single or few layers of 
gra phite0,i] allows to test the singular transport prop- 
erties predicted in early theoretical studies[l,|3,[l,[a| and 
experiments 0- The discovery of a substantial field ef- 
fect 8] and of ferromagnetic behavior allows to envis- 
age graphene as a reasonable replacement of nanotubes 
in electronic applications. Disorder plays a very impor- 
tant role in the electronic properties of low dimensional 
materials. 

Substitution of an hexagon by an n-sided ring in the 
hexagonal lattice without affecting the threefold coordi- 
nation of the carbon atoms leads to the warping of the 
graphene sheet. Rings with n > 6 (ri < 6), induce lo- 
cally positive (negative) curvature. Inclusion of an equal 
number of pentagons and heptagonal rings would keep 
the flatness of the sheet at large scales and produce a flat 
structure with curved portions that would be structurally 
stable and have distinct electronic properties. This de- 
fects give rise to long range modifications in the elec- 
tronic wave function that affect the electronics in a way 
different from that produced by vacancies or other impu- 
rities modelled by local potentials. Pentagon-heptagon 
pairs and Stone- Wales defects made of two adjacent hep- 
tagons and two pentagons form naturally in experiments 
of ion bombarded nanotubes as a mechanism to reduce 
the dangling bonds in large vacancies [l0| and have been 
observed to form in situ in single graphene layers by high- 
resolution transmission electron microscopy (TEM)[]J. 

We propose a method, based on a cosmological anal- 
ogy, to compute the electronic structure and transport 
properties of curved graphene sheets consisting of an ar- 
bitrary number of topological defects located at fixed po- 
sitions of the lattice. We see that, unlike vacancies and 
voids, the combination of positive and negative curva- 
ture gives rise to characteristic inhomogeneous patterns 
in the density of states that affect the transport proper- 
ties of the layers and can be observed in scanning tunnel 
(STM) and electron transmission spectroscopy (ETS). 
The present analysis can help in the experimental charac- 
terization of graphene samples and in the correct analysis 
of STM images. The results obtained can be related to 
recent Electrostatic Force Microscopy (EFM) measure- 
ments that indicate large potential differences between 



micrometer large regions on the surface of highly ori- 
ented graphite 1 121 and to the suppression of magnetorre- 
sistance in single layered graphene reported in [13] . We 
apply the formalism to present the corrections to the lo- 
cal density of states induced by pentagon-heptagon pairs 
and Stone- Wales defects in the average planar graphene 
sheet. 

Description of the model. The conduction band of 
graphene is well described by a tight binding model which 
includes the tt orbitals which are perpendicular to the 
plane at each C atom[l3, [l^. This model describes a 
semimetal, with zero density of states at the Fermi en- 
ergy, and where the Fermi surface is reduced to two in- 
equivalent K-points located at the corners of the hexago- 
nal Brillouin Zone. The low-energy excitations with mo- 
menta in the vicinity of any of the Fermi points K± have 
a linear dispersion and can be described by a continu- 
ous model which reduces to the Dirac equation in two 
dimensions [l6l [TtI . [isj . In the absence of interactions or 
disorder mixing the two Fermi points the electronic prop- 
erties of the system are well described by the effective 
low-energy Hamiltonian: 

Ho± = ihvY{(Txdx ± CTydy) , (1) 

where ax.y are the Pauli matrices, v-p = {3ta)/2, and 
a — l.AA is the distance between nearest carbon atoms. 
The components of the two-dimensional spinor: ^'i(r) — 
(93^(1"), 'PB(r))i correspond to the amplitude of the wave 
function in each of the two sublattices (A and B) which 
build up the honeycomb structure. Pentagons and hep- 
tagons can be viewed as disclinations in the graphene lat- 
tice, and, when circling one such defect, the two sublat- 
tices in the honeycomb structure as well as the two Fermi 
points are exchanged. The scheme to incorporate this 
change in a continuous description was discussed in refs. 
(TgI [l7j and [l^ . The process can be described by means 
of a non Abelian gauge field, which rotates the spinors in 
the flavor space of the Fermi points. The two spinors as- 
sociated to each Fermi point can be combined into a four 
component Dirac spinor 4'(r) — (^'+(r), \I'_(r)) which in 
the presence of a single disclination is an eigenstate of 
the Hamiltonian [2^ 

H = -iv^^.B + 57^7.^(0. (2) 
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Vp is the Fermi velocity, 7' are 4x4 matrices constructed 
from the Pauh matrices, j'' = ^ ^ I, the latin indices 
run over the two spatial dimensions and g is a coupling 
parameter related to the deficit angle of the defect. 

We use time-independent perturbation theory to cal- 
culate corrections to the self-energy in the weak coupling 
regime of the parameter g = , where the constant 
$ is the strength of the vortex: $ = ^ Adr related to 
the opening angle of the defect by $ = (1 — 5). L is the 
dimension of the sample. 

In a disordered system, in general, because of the pres- 
ence of an external potential or impurities, the space is 
inhomogeneous. The Green's function doesn't depend on 
the difference (r — r') and k is no longer a good quantum 
number. If the external potential or the effect of the im- 
purities are time-independent we have elastic scattering 
and the states k and k' have the same energy. We want 
to calculate the total density of states p{u!) of the system 
perturbed by the defect via the vector potential given in 
eq. (2). This density is the imaginary part of the Green's 
function integrated over all positions, in the limit r' r: 



J ImG{uj,r,r)dr. 



In terms of the Green's function in momentum represen- 
tation, p{lo) can be written as: 



dk 



dk' 



{2^? J (2^) 
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G{uj,k,k!)dv. 



The integration over r gives delta function 47r^(5(k — k'), 
and p{lo) then reads: 



p{uj) 



dk 

(2^ 



/mG(w, k,k). 



Generalization to negative curvature and an arbitrary 
number of defects. An alternative approach to the gauge 
theory of defects is to include the local curvature in- 
duced by an n-membered ring by coupling the Dirac 
equation to a curved space. This approach was intro- 
duced in [ll, [l3| to study the electronic spectrum of 
closed fuUerenes and has recently been used for fuUerenes 
of different shapes in 2l| . The holonomy discused in the 
gauge approach is included in the formalism by means of 
the spin connection [22|. In most of the previous works 
in this context, the graphene sheet was wrapped on a ge- 
ometrical surface that can be easily parametrized. The 
situation is more complicated if we try to describe the 
samples that are being obtained in the laboratory that 
are flat sheets with local corrugation [Tsj. 

In this context one can see that the substitution of an 
hexagon by a polygon of n < 6 sides gives rise to a conical 
singularity with deficit angle (27r/6)(6 — n). This kind of 
singularities have been studied in cosmology as they are 
produced by cosmic strings, a type of topological defect 
that arises when a U(l) gauge symmetry is spontaneously 
broken [2^ . We can obtain the correction to the density of 
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FIG. 1: Electronic density around a conical defect. 



states induced by a set of defects with arbitrary opening 
angle by coupling the Dirac equation to a curved space 
with an appropriate metric as described in ref. [l^l . The 
metric of a two dimensional space in presence of a single 
cosmic string in polar coordinates is: 



ds' 



-dt^ 



dr^ + c^r^de^, 



(3) 



where the parameter c is a constant related to the deficit 
angle by c = 1 — 6. 

The dynamics of a massless Dirac spinor in a curved 
spacetime is governed by the Dirac equation: 



(4) 



The difference with the flat space lies in the deflnition of 
the 7 matrices that satisfy generalized anticommutation 
relations 



{7^7^} = 25'^^ 



(5) 



where g^^ is given by ([3]) , and in the covariant derivative 
operator, deflned as 



(6) 



where F^^ is the spin connection of the spinor field that 
can be calculated using the tetrad formalism [25|. Fig. 
[1] shows the solution of the Dirac equation (jH) in the 
presence of a single defect with a positive deficit angle 
(positive curvature). The electronic density is strongly 
peaked at the position of the defect suggesting a bound 
state but the behavior at large distances is a power law 
with an angle-dependent exponent less than two which 
corresponds to a non normalizable wave function. This 
behavior is similar to the one found in the case of a single 
vacancy [26] and suggests that a system with a number 
of this defects with overlapping wave functions will be 
metallic. 

The case of a single cosmic string which represents a 
deficit angle in the space can be generalized to describe 
seven membered rings representing an angle surplus by 
considering a value for the deficit angle c greater than 1. 
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This situation is non-physical from a general relativity 
viewpoint as it would correspond to a string with nega- 
tive mass density but it makes perfect sense in our case. 
The scenario can also be generalized to describe an arbi- 
trary number of pentagons and heptagons by using the 
following metric: 

where A(r) = Y.f=i log(ri) and n = [{x - a^)^ + (y - 
bi)^]^^'^- This metric describes the space-time around N 
parallel cosmic strings, located at the points (a^, bi). The 
parameters are related to the angle defect or surplus 
by the relationship Ci = 1 — 4/ii in such manner that if 
Ci < 1(> 1) then > 0(< 0). 

From equation Q we can write down the Dirac equa- 
tion for the electron propagator, S'i?(a;, a;'): 

ij^{r)idf, - r^)SF{x, x') = S^{x - x'), (8) 

where x = {t,r). The local density of states N{LU,r) 
is obtained from ^ by Fourier transforming the time 
component and taking the limit r r': 

N{uj, r) = ImTrSpiuJ, r, r). (9) 

We solve eq. ([9]) considering the curvature induced by 
the defects as a perturbation of the flat graphene layer. 
The details of the calculation will be given elsewere[20|. 
Here we will show the results obtained. 

Results. Fig. [5] shows the correction to the local den- 
sity of states at energy w — 2.1eV and for a large region 
of the graphene plane with a pentagon (p)-heptagon (h) 
pair located out of plane at p = [—0.3, —5], h = [0.3, —5]. 
The LDoS is normalized with the clean DoS of graphene 
at the given energy. The color code is indicated in the fig- 
ure: green stands for the DOS of perfect graphene at the 
given energy and red (blue) indicates an accumulation 
(depletion) of the density in the area. The correction 
obtained is of the order of a few percent. We can see 
that pentagonal (heptagonal) rings enhance (deppress) 
the electron density. A similar result has been obtained 
in [27} with numerical simulations. It is to note that a 
somehow contradictory result was obtained in [28| where 
they studied the electrostatics of a graphene plane with 
defects. They found that disclinations corresponding to 
rings with more (less) than six carbon atoms function 
as attractors (repellent) to point charges. It is obvious 
that this issue needs further investigation. The dipo- 
lar character of the defect is clear from the figure and 
makes this type of defects unique. The effects produced 
will extend spatially much beyond the scale of the de- 
fect what can be related with recent Electrostatic Force 
Microscopy (EFM) measurements that indicate large po- 
tential differences between micrometer large regions on 
the surface of highly oriented graphite that do not show 
inhomogeneities in STM 12]. The interference patterns 
in the DOS depend on the position and on the relative 
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FIG. 2: (Color online) Image of the correction to the lo- 
cal density of states in a large portion of the plane with 
a heptagon-pentagon pair located out of plane at position 
[0.3, -5], [-0.3, -5] . 




FIG. 3: (Color online) Correction to the local density of states 
around a Stone- Wales defect located out of the image in the 
upper right corner. 



orientation of the dipoles. Fig. [3] shows the patterns pro- 
duced by a Stone- Wales defect (two pairs of pentagons 
(pi) heptagons (hi)) located out of plane at positions 
Pi = [5,5.6],p2 = [6.5,6], hi = [6,5.5],/i2 = [6,6.9] for 
a frequency w = 2.1eV. 

A random distribution of defects. A random distribu- 
tion of topological defects can be described in the contin- 
uous model by a (non abelian) random gauge field. The 
disorder is defined by a single dimensionless quantity. A, 
which is proportional to the average fluctuations of the 
field: 

(A(r)A(r')) - A(52(r-r') (10) 

The statistical properties of the gauge field induced by 
topological defects were analyzed in ref. [l^ following a 
renormalization group scheme and a complete phase di- 
agram as a function of disorder and Coulomb interaction 
was obtained in [2ff]. 

In ref. [19] it was shown that when pentagons and hep- 
tagons are bound in dislocations with average distance 6, 



4 



the vector field of a vortex-antivortex dipole decays as 
and the behavior of A was found to be: 

A cx ^Inb^ (11) 

where n is the density of dislocations. 

Assuming that random fields induced by topological 
defects have the same statistical properties to those with 
gaussian disorder with the same value of A, a renormal- 
ization group (RG) analysis using the replica trick gives 
rise to an effective interaction between fermion fields in 
different replicas. The resulting self-energy is logarith- 
mically divergent what can be interpreted as a renormal- 
ization of the density of statesfioj. The results obtained 
in this paper show a qualitative picture of the corrected 
DoS. The influence of these defects on the localization in 
graphene has been analyzed recently in [3d| . 

The type of disorder analyzed in this article belongs 
to the random gauge field in the classification of [23|. 
There it was shown that the combination of this type of 
disorder with the unscreened Coulomb interaction give 
rise to a line of infrared stable fixed points. Topological 
defects give rise also to long range correlated disorder 
whose properties will be analyzed else where [20l|. 

Conclusions and open problems. We have presented in 
this article results on the electronic structure of a curved 
graphene sheet with pentagon-heptagon pairs that we ex- 
pect to be present in the samples that are subject of 
present investigation. The formalism can also apply to 
clean samples of graphite showing strong anisotropy and 



to curved samples of graphene with an arbitrary num- 
ber of defects with any sign of the curvature. The model 
can also be used to study transport properties of the 
curved graphene sheets and of nanographite and results 
will be presented in The main achievement of this 
paper is the proposal of a formalism that can describe 
the corrugations observed in the experiments and pre- 
dicts charge inhomogeneities that can be used to charac- 
terize the samples. We have computed the local density 
of states with the geometrical formalism found that it is 
enhanced around defects which induce positive curvature 
in the lattice while the charge is "repelled" from regions 
with negative curvature. Heptagon-pentagon pairs that 
keep the graphene sheet flat in the long range behave 
as dipoles and give rise to characteristic modulations of 
the DOS that can be observed by STM and can explain 
recent experiments. The magnitude of the oscillations in- 
creases with the frequency and the characteristic patterns 
could be used to characterize the graphene samples. The 
features predicted in this work should also be observable 
in other layered materials with similar structure as boron 
nitride [3l|. The present analysis can help to clarify the 
issue of the analysis and interpretation of STM images. 
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